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Abstract We study one of the simplest covariant modified-gravity models based on 
the Dvali-Gabadadze-Porrati (DGP) brane cosmology, a self-accelerating universe. In 
this model gravitational leakage into extra dimensions is responsible of late-time accel- 
eration. We mainly focus on the effects of the model parameters on the geometry and 
the age of universe. Also we investigate the evolution of matter density perturbations 
in the modified gravity model, and obtain an analytical expression for the growth in- 
dex, /. We show that increasing Q rc leads to less growth of the density contrast 8, and 
also decreases the growth index. We give a fitting formula for the growth index at the 
present time and indicate that dominant term in this expression verifies the well-known 
approximation relation / ~ Q^i- As the observational test, the new Supernova Type 
la (SNIa) Gold sample and Supernova Legacy Survey (SNLS) data, size of baryonic 
acoustic peak from Sloan Digital Sky Survey (SDSS), the position of the acoustic peak 
from the CMB observations and the Cluster Baryon Gas Mass Fraction (gas) are used 
to constrain the parameters of the DGP model. We also combine previous results with 
large scale structure formation (LSS) from the 2dFGRS survey. Finally to check the 
consistency of the DGP model, we compare the age of old cosmological objects with 
age of universe in this model. 

PACS numbers: 05.10.-a ,05.10.Gg, 05.40.-a, 98.80. Es, 98.70.Vc 



1 Introduction 

Recent Observations of type la supernova (SNIa) provides the main evidence for accel- 
erating expansion of the Universe [HE]. Analysis of SNIa and the Cosmic Microwave 
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Background radiation (CMB) observations indicates that about 70% of the total energy 
of the Universe is made by the dark energy and the rest of it is the dark matter with a 
few percent of Baryonic matter 3,45 . The " cosmological constant" is a possible ex- 
planation for the acceleration of the universe [6]. This term in Einstein field equations 
can be regarded as a fluid with the equation of state of w = — 1. However, there are 
two problems with the cosmological constant, namely the fine-tuning and the cosmic 
coincidence. In the framework of quantum field theory, the vacuum expectation value 
is 123 order of magnitude larger than the observed value of 10™ 47 GeV 4 . The absence 
of a fundamental mechanism which sets the cosmological constant to zero or to a very 
small value is the cosmological constant problem. The second problem known as the 
cosmic coincidence, states that why are the energy densities of dark energy and dark 
matter nearly equal today? 

There are various solutions for this problem as the decays cosmological constant 
models. A non-dissipative minimally coupled scalar field, the so-called Quintessence 
model can play the role of this time varying cosmological constant [7][8][9]. The ratio 
of energy density of this field to the matter density in this model increases by the 
expansion of the universe and after a while dark energy becomes the dominant term 
of the energy-momentum tensor. One of the features of this model is the variation of 
equation of state during the expansion of the universe. Various Quintessence models like 
k-essence [10], tachyonic matter [11], Phantom [121113] and Chaplygin gas [14] provide 
various equations of states for the dark energy [T3irT4irT5T[T6Tl 1 7111 8 l[T9l[20l[2T| . 

Another approach dealing with this problem is using the modified gravity by chang- 
ing the Einstein-Hilbert action. Some of models as 1/7? and logarithmic models provide 
an acceleration for the universe at the present time [22) . In addition to the phenomeno- 
logical modification of action, the brane cosmology also implies modification for the 
general relativity on a brane embedded in an extra dimension space. Some brane world 
models which produce the late time acceleration have been tested using many obser- 
vational experiments such as local gravity 23; 24 ,25 , Supernova Type la 22,26,27 
28,29,30,31,32 , angular size of compact ratio sources [33], the age measurements of 
high redshift objects [34], the optical gravitational lensing surveys [35], the large scale 
structures [36], and the X-ray gas mass fraction in galaxy clusters [37U38] . 

In some recent papers 31,39,40 observational constraints have been obtained 
through the old data of Supernova Gold sample and its combination with CMB shift 
parameter and Baryon acoustic oscillation. Recently Guo et al. have put constraints 
on this model using recent SNIa data and Baryon acoustic oscillation [26]. Song et al., 
Sawichi and Carroll have separately investigated the effect of DGP on the integrated 
Sachs- Wolfe and tested the validity of modified linear growth factor in the sub-horizon 
scale [411142] , In ref. [41] linear growth of density contrast in this model has been re- 
viewed but some of the interesting quantity as the growth index and its behavior versus 
redshift and its dependency to the model parameter is missed. 

In this paper we examine the effects of DGP model on the geometrical parameters 
of the universe. On the other hand we use the observational results related to the back- 
ground evolution. Since the structure formation in DGP is currently well understood 
on scales between a few percent of the Hubble scale and the scale radius of a typical 
dark matter halo [44], we combine those results with the linear structure formation 
of large scale in the universe. Meanwhile we concentrate our attention to the effect 
of Q Tc and Q m as free parameters of the model on the density contrast and growth 
index evolution. We extend the simplest growth index analytic formula given for the 
flat ylCDM [45] in the underlying modified gravity model. We organize this paper as 
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follows: In Sec[TJ we introduce DGP model as a self-accelerating cosmology. Its free 
parameters and modified Friedman equation which governs on the background dynam- 
ics of the universe are also investigated. In Sec[2]we study the effect of this model on 
the comoving distance, comoving volume element, the variation of angular size by the 
redshift [46], In Sec. [3] we put some constraints on the parameters of model by using 
the background evolution, such as new Gold sample and Legacy Survey of Supernova 
Type la data |47| . the position of the observed acoustic angular scale on the last scat- 
tering surface, CMB shift parameter, the baryonic oscillation length scale and baryon 
gas mass fraction for the range of redshift, z < 1.0. We study the linear structure 
formation in this model and compare the growth index with the observations from the 
2— degree Field Galaxy Redshift Survey (2dFGRS) data in Sec. [4] We also compare the 
age of the universe in this model with the age of old cosmological structures in Sec. [5] 
Sec|6] contains summary and conclusion of this work. 



2. DGP MODIFIED GRAVITY 

One of the simplest covariant modified-gravity models is based on the Dvali-Gabadadze- 
Porrati (DGP) brane- world model, as generalized to cosmology by Deffayet |48| . (It is 
worth noting that the original DGP model with a Minkowski brane was not introduced 
to explain acceleration - the generalization to a Friedman brane was subsequently 
found to be self-accelerating.) In this model, gravity leaks off the 4-dimensional brane 
universe into the 5-dimensional bulk spacetime at large scales. Ordinary matter is 
considered to be localized on the brane while gravity can propagate in the bulk. At 
small scales, gravity is effectively bound to the brane and 4D gravity is recovered to a 
good approximation. The action for the five-dimensional theory is: 

S = / d°X ^-3(5)^(5) + / d 4 2^-5(4)^(4) + ^matter (1) 

zk 5 j zk 4 j 

where the subscripts 4 and 5 denote the quantities on the brane and in the bulk, 
respectively, k^(k^) is the inverse of four(five)-dimensional reduced Planck mass, and 
Smatter is the action for matter on the brane. The solution of DGP action in FRW 
metric provides a self-accelerating universe for the expanding phase of the universe. 
This model can be an alternative to the cosmological constant for describing the present 
acceleration of the universe. However, DGP model suffers from the ghost instability 
that was shown in 49,50 though the boundary effective action formalism. On the other 
hand the existence of ghost is confirmed by explicit calculation of the spectrum of linear 
perturbations in the five-dimensional framework 51 . A solution for this problem is so- 
called cascading DGP model in which the unlike previous attempts, it is free of ghost 
instabilities. In this model the 4D propagator is regulated by embedding the 3-brane 
within a 4-brane with their own gravity terms induced by a flat 6D bulk [52] , 

Coming back to the action (fT|, at moderate scales the induced gravity term is 
responsible for the recovery of 4-Dimensional Einstein gravity. The transition from 4D 
to 5D behavior is governed by a crossover scale r c : 

In the weak-field gravitational field, potential behaves as r _1 for r <C r c and as r~ 2 
for r 3> r c . At large scale gravity is five dimensional. The energy conservation equation 
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remains the same as in general relativity, but the Friedman equation is modified: 

p + 3H(p + p) = 0, (3) 

(4) 



rr 2 K . 1 / rr9 K 
H + ± —\ H 2 + 

a £ r c V a 



8nG 



Two given sign in equation @ correspond to the two branches of the cosmological 
evolution. The upper sing shows a de Sitter expansion of the universe, while the lower 
sign corresponds to the self-accelerating solution without the cosmological constant. 
So we infer the cosmological effects of the second branch of DGP model. Equations © 
and Q imply (for the CDM case p = 0) 



H 



K 



= -47TG/5 



yji + 32TvGr 2 p/3 



(5) 



Equation Q shows that at early times, when H 2 + K/a 2 ^> r^ 2 , the general relativistic 
Friedman equation is recovered. By contrast, at late times in a CDM universe, with 
' 1 0, we have 



p oc a 



H 



Hoc = — 



(6) 



Gravity leakage at late times initiates acceleration not due to any negative pressure 
field, but due to the weakening of gravity on the brane. Since Hq > Hoc, in order to 
achieve self-acceleration at late times, we require 



and this is confirmed by fitting observations as discussed below. 

In dimensionless form, the modified Friedman equation (j4| is given by 



(J) 



H 2 



+ n K (i + zf 



'Or 



(8) 



where 



f^fC — 1 — ^rn — 2 a/ £2 r 



(y/0i~+ y/Qr c + ft m ) 



1 — fitot 
1 



(9) 
(10) 



From equation ([5]), the dimensionless acceleration is 
1= jp^ = (y«m(l + ^) 3 + «r c + 



'Or 



+ 



2Q ra - OmQ- + zf 

2yjn m (i + z) 3 + n r 



(ii) 
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so that the redshift at which acceleration era is started is given by [37] 



Z„=Q = 2 



Also equation <(9j shows in the flat model: 

n rc = j(i~n m f (13) 

For the transition from negative to the positive acceleration at the present-time, equa- 
tion (|12l) implies : 

Or. = ^ (14) 

Acceleration parameter for this model in terms of scale factor is shown in Figure [1] 
Increasing the value of Q Tc causes that universe entered in the acceleration phase at 
the earlier times. The lower panel of Figure [T] shows acceleration parameter for the flat 
ylCDM model, obviously fi Tc has the same role as cosmological constant. 

The modified Friedman equation in DGP may be reinterpreted from a standard 
viewpoint. We define the effective dark energy density p e ff = 3H /8irGr c . Then the 
effective dark energy equation of state W e ff = Peff/ Peff is given by p e ff + 377(1 + 
w e ff)p e ff = 0. Thus Peff and w e jf give a standard general relativistic interpretation 
of DGP expansion history, i.e., they describe the equivalent general relativity dark 
energy model. For the flat case, fix = 0, we find 



f . n m -l-J(l- J? m ) 2 + 4J? m (l + 2 )3 

W e tf(Z)= ; , (15) 

which implies 

^//(°) = -irk- (16) 

The DGP and ylCDM models have the same number of parameters, with r c sub- 
stituting A, therefore DGP model gives a very useful framework for comparing the 
ylCDM general relativistic cosmology to a modified gravity alternative. Now an inter- 
esting question that arises is: "can DGP model predict dynamics of universe?" or in 
another word, "what values of the model parameter to be consistent with observational 
tests?" 

In the forthcoming sections we will study the observational constraints on the 
model. 



2 The effect of DGP model on the geometrical parameters of universe 

The cosmological observations are mainly affected by the background dynamics of 
universe. In this part we study the sensitivity of the geometrical parameters on the 
parameters of DGP model. 
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Fig. 1 Upper panel shows acceleration parameter (q = a/H^a) in the DGP model as a 
function of scale factor for various values of Q ra . Lower panel corresponds to the same function 
for the flat ylCDM. We chose the flat universe. 



2.1 comoving distance 

The radial comoving distance is one of the basic parameters in cosmology. For an 
object with the redshift of z, using the null geodesies in the FRW metric, the comoving 
distance is obtained as: 

r(z-,n m ,n rc ) = )=t(^\o^\ f d f ) , 

HoV\^\ V Jo H (z')/ H 0j 

(17) 



£2 r =0.05 
O r c =0.10 
=015 

£V =0.20 

c 





Observational Constraints with recent Data on the DGP Modified Gravity 



7 




Fig. 2 Comoving distance, r(z; Q m , Q Tc ) (in unit of c/ Hq) as a function of redshift for various 
values of i?,- e . We fixed Qjf = 0.0. 



where 

J-(x) = (x, sin x, sinh x) for K — (0, 1, — 1) . (18) 

and H(z; n m , flr c ) is given by equation ([8]). By numerical integration of equation 
(|17[) . the comoving distance in terms of redshift for different values of f2 Tc is shown 
in Figure [5] Increasing the O rc results in a longer comoving distance. According to 
this behavior, by tuning the value of Q Tc we may explain the Supernova Type la 
observations. 



2.2 Angular Size 

The apparent angular size of an object located at the cosmological distance is another 
important parameter that can be affected by the cosmological model during the history 
of the universe. An object with the physical size of D is related to the apparent angular 
size of 6 by: 

D = d A 6 (19) 

where d A = r(z; Q m , O rc )/(l + z) is the angular diameter distance. The main applica- 
tions of equation (I19|l is on the measurement of the apparent angular size of acoustic 
peak on CMB and baryonic acoustic peak at the high and low redshifts, respectively. 
By measuring the angular size of an object in different redshifts (the so-called Alcock- 
Paczynski test) it is possible to probe the validity of modified gravity models [46]. The 
variation of apparent angular size AO in terms of Az is given by: 

^ = H(z;n m ,n rc )r(z;n m ,n rc ) (20) 



s 
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Fig. 3 Alcock-Paczynski test comparing Az/AO as a function of redshift for four different 
flr c normalized to the case with Q m = and fl Ta = 0.25 (fiat universe fife = 0). It must be 
pointed out that for other values of fl Ta we also assumed fix = 0.0. 



Figure [3] shows Az/A6 in terms of redshift, normalized to the case with Q m = 0.0 
and flat universe fix — 0.0. The advantage of Alcock-Paczynski test is that it is 
independent of standard candles and knowing a standard ruler such as the size of 
baryonic acoustic peak one can use it to constrain the modified gravity model. 



2.3 Comoving Volume Element 

The comoving volume element is another geometrical parameter which is used in 
number-count tests such as lensed quasars, galaxies, or clusters of galaxies. The co- 
moving volume element in terms of comoving distance and Hubble parameter is given 
by: 

/0; Q m , Q Tc ) ee = r 2 (z; fl m , f2 Tc ) / H(z; fl m , f?r e ). (21) 

According to Figure 21 the comoving volume element becomes large for larger value of 
fl Tc in the flat universe. 



3 Observational Constraints From the Background Evolution 

In this section we compare the SNIa Gold sample released from more recent observa- 
tions 53 which have lower systematic errors than the former Gold sample data set. 
This new catalog contains 156 Supernova type la. On the other hand we also take into 
account 116 supernova Legacy Survey data which seems to be more consistent with 
WMAP observation as another SNIa observations to examine DGP model. To make 
the acceptance interval of model free parameters more confined, we use the location 
of acoustic peak of temperature fluctuations from WMAP observation, the location of 
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Fig. 4 The comoving volume element in terms of redshift for various fl Tc exponent. Increasing 
Q Tc shifts the position of the maximum value of volume element to lower redshifts. We fixed 
i? K =0.0 



baryonic acoustic oscillation peak from the SDSS and baryon gas mass fraction in the 
cluster for 26 samples at redshifts less than 1 [54] . The Supernova Type la experiments 
provided the main evidence of the existence of dark energy. Since 1995 two teams of 
the High-Z Supernova Search and the Supernova Cosmology Project have discovered 
several type la supernovas at the high redshifts [161155] . Recently Riess et al.(2004) 
announced the discovery of 16 type la supernova with the Hubble Space Telescope. 
This new sample includes 6 of the 7 most distant (z > 1.25) type la supernovas. They 
determined the luminosity distance of these supernovas and with the previously re- 
ported algorithms, obtained a uniform 156 Gold sample of type la supernovas as a new 
data set with lower systematic errors than former Gold sample data [47|I56| I57]. At the 
beginning we compare the predictions of the DGP model with the recent SNIa Gold 
sample [53] . The observations of supernova measure essentially the apparent magni- 
tude m including reddening, K correction, etc, which are related to the (dimensionless) 
luminosity distance, D^, of an object at redshift z through: 

m = M + 51ogD L (z; f2 m , ftr c ), (22) 

where 

Di(=:Bm - aj ^K^Iwi)' <23) 

Also 

M = M + 51og( -g^g ) +25. (24) 
y 1 Mpc J 

where M is the absolute magnitude. The distance modulus, /i, is defined as: 
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(i = m — M = 5 log D^(z; O m , O, 

c/H 
1 Mpc 



+ 51 °g( l^TT- ) + 25 > ( 25 ) 



H = 5]ogD L (z;f2 m ,n r<s ) + M (26) 

In order to compare the theoretical results with the observational data, we must 
compute the distance modulus, as given by equation (|25|l . For this purpose,the first 
step is to compute the quality of the fitting through the least squared fitting quantity 
X defined by: 

2/-lff o O \ _ [Mofcs( 2 i) - Mth(z»; ^m,^r c ,M)] 2 
X (iM ,Si m ,!i r ) — > 5 , 

i 

(27) 

where <Tj is the observational uncertainty in the distance modulus. To constrain the 
parameters of model, we use the Likelihood statistical analysis: 

C(M,n m ,n rc ) =Me- x2{M ' nm ^ c )/2 (28) 

where M is a normalization factor. The parameter M is a nuisance parameter and 
should be marginalized (integrated out) leading to a new \ 2 defined as: 

r+oo 

X 2 = -2hi/ C{M,n m ,n re )dM (29) 



Using equations ([27) . ((28)) and ((29)) . we find 



X (tt m , n r J =X (M = 0, J?m, ftrj ^ — 

+ ln(C/27r) (30) 

where 

B(a m , n rc ) = ^^)-^(^^-^,M = o)] | (31) 



C = £4 (32) 



Equivalent to marginalization is the minimization with respect to M. One can show 
that x" 2 can be expanded in M as [61] : 

X 2 {O m , Tc ) = X 2 (M = 0, Q m , (2r c ) - 2MB + M 2 C* (33) 
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Fig. 5 Distance modulus of the SNIa new Gold sample in terms of redshift. Solid line shows 
the best fit values with the corresponding parameters of h = 0.64, Q m = 0.36^Q Qg, fl r< . = 

0.23l°'ol in lc level of confidence with xLn/ W <i.»./ = 091 for DGP model. 




Fig. 6 Distance modulus of the SNLS supernova data in terms of redshift. Solid line shows 
the best fit values with the corresponding parameters of h = 0.70, Q m = 0.13^q'q|, Q rc = 
0-14+003 in la level of confidence with X 2 min /N d . . f = 0.85 for DGP model. 

which has a minimum for M = B /C: 

XSNI^m, Or c ) = X\M = 0, Q m , Or a ) - B ^ m ^ r ^ (34) 

Using equation (|34|) we can find the best fit values of model parameters as the values 
that minimize x 2 (^m, O r< J). For the Likelihood analysis we use some weak priors for the 
model parameters indicated in Table[T] The best values for the parameters of the model 
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Table 1 Priors on the parameter space, used in the likelihood analysis. 



Parameter 


Prior 




Otot 




Free 




0.00 - 1.00 


Top hat 




0.00 - 1.00 


Top hat 




0.020 ± 0.005 


Top hat (BBN) 


h 







are: Q m = 0.36±°;°£, Q Tc = 0.23±°;°4 and Q K = -0.56^ with xLiJNd.cf = 0.91 
at la level of confidence. The corresponding value for the Hubble parameter at the 
minimized \ 2 ls h = 0.64 and since we have already marginalized over this parameter 
we do not assign an error bar for it. The best fit values for the parameters of model by 
using SNLS supernova data are !2 m = 0.13+°°^, ft Ta = 0.14^03 and Q K = 0.20±°^ 
with Xmin/Nd.o.f = 0-85 at la level of confidence. The value of Hubble parameter at 
the minimum value of \ 2 is h = 0.70. Obviously our results are different from what 
reported in |26| and |31| . In the first reference they report the following values for the 
model parameters: Q m = 0.34^q'q 8 and Q Ta = 0.24^q'q^ using old SNIa Gold sample 
while in the second reference Marteens et. al. reported: (2 m = 0.270 and !2 rc = 0.125. 
We point out that the constraint by SNIa is very sensitive to the various catalogs of 
supernova data set 62 . Table [2] indicates the results from observational constraints on 
the free parameters. Figures [S] and [5] show the comparison of the theoretical prediction 
of distance modulus by using the best fit values of model parameters and observational 
values from new Gold sample and SNLS supernova, respectively. For the age consistency 
test we substitute the parameters of model from the SNIa new Gold sample and SNLS 
fitting in equation (|66|) (see Sec.[S]for more details) and obtain the age of universe about 
13.78^Q'gg Gry and 14. 96^ 1 ' 43 Gry, respectively. They give a universe older than what 
is expected from the old stars. 

The other constrain results from the CMB acoustic peak observations. Before last 
scattering, the photons and baryons are tightly coupled by Compton scattering and 
behave as a fluid. The oscillations of this fluid, occurring as a result of the balance 
between the gravitational interactions and the photon pressure, lead to the familiar 
spectrum of peaks and troughs in the averaged temperature anisotropy spectrum which 
we measure today. The odd and even peaks correspond to maximum compression of 
the fluid and to rarefaction, respectively |64] , In an idealized model of the fluid, there 
is an analytic relation for the location of the m-th peak: l m ~ ml a [651166] where I a 
is the acoustic scale which may be calculated analytically and depends on both pre- 
and post-recombination physics as well as the geometry of the universe. The acoustic 
scale corresponds to the Jeans length of photon-baryon structures at the last scattering 
surface some ~ 379 Kyr after the Big Bang [5]. The apparent angular size of acoustic 
peak can be obtained by dividing the comoving size of sound horizon at the decoupling 
epoch r s (zdec) by the comoving distance of observer to the last scattering surface 

r(zdec)- 

9 A = f= (35) 
Ia r(z dec ) 



The size of sound horizon at the numerator of equation (|35[) corresponds to the distance 
that a perturbation of pressure can travel from the beginning of universe up to the last 
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scattering surface and is given by: 

where v s (z)~ 2 = 3 + 9/4 x pb(z) / p ra d( z ) is the sound velocity in the unit of speed of 
light from the big bang up to the last scattering surface [181165] and the redshift of the 
last scattering surface, z dec , is given by [53] : 

z dec = 1048 [1 + 0.00124Kr a738 ] [1 + gi (io m ) 92 ] , 
gi = 0.0783Kr ' 238 [l + 39.5K) ' 763 ] _1 , 

32 = 0.560 [l + 21.1( W6 ) 1 - 81 ]"\ (37) 

where tum = Omh 2 , = Q^h 2, and p ra( i is the radiation density. is relative baryonic 
density to the critical density at the present time. Changing the parameters of the 
model can change the size of apparent acoustic peak and subsequently the position of 
I A = tt/^a in the power spectrum of temperature fluctuations at the last scattering 
surface. The simple relation l m m ml a however does not hold very well for the peaks 
although it is better for higher peaks [66II67| . Driving effects from the decay of the 
gravitational potential as well as contributions from the Doppler shift of the oscillating 
fluid introduce a shift in the spectrum. A good parametrization for the location of the 
peaks and troughs is given by [66H67] 

lm = lA(m- 4>m) (38) 

where <f> m is phase shift determined predominantly by pre-recombination physics, and 
are independent of the geometry of the Universe. The location of acoustic peaks can 
be determined in model by equation (|38[l with <f> m {u) m , Doran et. al. [67], recently 
have shown that the first and third phase shifts are approximately model independent. 
The values of these shift parameters have been reported as: <f>i(ui m ,Ub) — 0.27 and 
4> 3 {Lo m ,Ld b ) ~ 0.341 [661167] . According to the WMAP observations: h = 220.1 ± 0.8 
and Z3 = 809 ± 7, so the corresponding observational values of read as: 

l°A S \h = (1 h (j)i) = 299.45 ± 2.67 (39) 

^ bS l' 3 = (3 _^ 3) = 304 - 24 ± 2 - 63 ( 4 °) 
their Likelihood statistics are as follows: 



Xh = 2 ( 41 ) 

a l 

and 



r/obsi 7th 1 

l l A \h - A \h 



xi 3 = -2 ( 42 ) 

O o 



r/obsi 7th 1 
[ l A Us _^A \h 

4 



1 2 
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because of weak dependency of phase shift to the cosmological model one can use 
another model independent parameter which is so-called shift parameter 1Z: 

,flat 

TZ(x^—, (43) 
n 

where l^ lat corresponds to the flat pure-CDM model with £2 m =1.0 and the same io m 
and uif, as the original model. It is easily shown that shift parameter is as follows [68] : 

n = ^ D L^dec,n m n rc ) 

(1 + Zdec) 

The observational results of CMB experiments correspond to a shift parameter of 
U = 1.716 ± 0.062 (given by WMAP, CBI, ACBAR) [5]|69]. One of the advantages of 
using the parameter 1Z is that it is independent of Hubble constant. In order to put 
constraint on the model from CMB, we compare the observed shift parameter with 
that of model using likelihood statistic as [68] : 

C ~ e" x ™ B/2 (45) 

where 

°"CMB 

where TZ^ and TZ \^ S axe determined using equation (]44p and given by observation, 
respectively. Figures [7] shows Q Tc , and Oj\ in DGP model and the /1CDM as a function 
of f} m for a given Decreasing both f2 Ta and fi\ lead an increasing in the value of 
present matter density. 

2 2 -4-0 03 

Using XSNIa + Xu analysis, we find the best fit values as: Q m = 0.15_q' O3 , flr c = 
°- 08 io'o2 imPty n K = 0.42+° ^ with X 2 m in/ N d.o.f = 0.96 for New catalog of Gold sam- 
ple. For SNLS SNIa combined with the position of first peak we find: O m = O.H^q'q,, 

S7r c = O.lliSloi ^P 1 ^ Q K = 0.36i°;o5 with xLin/ N d.o.f = 0.85 at la level of confi- 
dence. Table [3] gives the best values of parameters using the location of first and third 
peaks of CMB power spectrum and other observational tests. We see that Q m and fl rc 
are very sensitive to the peaks position of power spectrum of temperature fluctuations 
at the last scattering surface. Since the phase transitions of peak position are weakly 
model dependent we also apply the shift parameter of CMB to extract best values for 
model parameters. According to the XsNIa + Xcmb statistic we get: J? TO = 0.23^q'q3, 
f2r c = 0.14+°;°J imply Q K = 0-03+g;g| with Xmin/ N d.o.f = 0.93 for new catalog 
of Gold sample. For SNLS SNIa combined with the position of first peak we find: 
n m = 0.17™, O rc = 0.16lg;gJ imply Q K = 0.05±°;°| with X 2 mln /Nd.o.f = 0.84 
at la level of confidence. The corresponding age of the universe are 15.23;°;^ and 
14.48^g'^g, respectively. These values are slightly different than that of reported in the 
previous analysis [26,31 . 

The recently detected size of baryonic peak in the SDSS is the third observational 
data for our analysis. The correlation function of 46,748 Luminous Red Galaxies (LRG) 
from the SDSS shows a well detected baryonic peak around 100/i -1 Mpc. This peak 
was identified with the expanding spherical wave of baryonic perturbations originating 
from acoustic oscillations at recombination. This peak has an excellent match to the 
predicted shape and the location of the imprint of the recombination-epoch acoustic 
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£2 m 

Fig. 7 Constant acoustic angular scale in the joint space of Q m and Q Ta (upper panel). Lower 
panel shows dependence of acoustic angular scale on the i? m and cosmological constant. 

oscillation on the low-redshift clustering matter [70]. Recently Linder has shown in de- 
tail some systematic uncertainties for baryon acoustic oscillation [71]. Nonlinear mode 
coupling which is related to this fact that the baryon acoustic oscillation is mostly 
contributed by linear scale, but the influence of non-linear collapsing has quite broad 
kernel. In other words, one might say that baryon acoustic oscillation are 90 — 99% 
linear in comparison to the CMB which is 99.99% linear, so this difference may affect 
on various models in different way. Careful works to constrain on the free parame- 
ters of underlying model needs to be carried out to determine the effect of nonlinear 
mode coupling in the results of constraint by SDSS observation. Nevertheless, roughly 
speaking regards the acceptance intervals for free parameter cover the real intervals 
determined by assuming nonlinearity mode for SDSS observation [72] . 
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Table 2 The best values for the parameters of DGP self accelerating model with the corre- 
sponding age for the universe from fitting with SNIa from new Gold sample and SNLS data, 
SNIa+CMB+gas, SNIa+CMB+gas+SDSS and SNIa+CMB+gas+SDSS+LSS experiments at 
one and two a confidence level. 



Observation 






n K 


Age (Gyr) 


SNIa(ncw Gold) 


o.36t° :E 
0.3618:12 


23+ - 04 
u - zo -0.04 

n 9 o+0.09 

u - zo -o.io 


-o 56+ - 20 

U.OD_ 20 

-n 56+ - 43 

U.OD_q 4g 


13.78l°-j« 


SNIa(ncw Gold) 
+CMB+gas 


o 2S+ - 04 

23+ 08 


14+0- 02 
u - ±y: -0.02 

14+0-03 
u - ±y: -0.03 


u - u -o.io 
04+ ' 18 

u ' u -0.16 


15.23l°;® 4 


SNIa(new Gold) 
+CMB+SDSS 
+gas 


30+ - 01 

u - ou -0.05 

30+ - 04 
u - ,3u -0.06 


12+0' 02 
u - lz -o.oi 

12+ ' 03 
u - ±z -0.03 


oi+o- 09 

u - ul _0.09 
01+0-15 


14 52+ - 15 


SNIa(ncw Gold) 
+CMB+SDSS 
+LSS+gas 


n 9 o+0.03 
n 9 o+0.05 


n i3+o- 01 
u - 1,3 -o.oi 

n i3+o- 02 

U - ±J -0.03 


-0 002+ ' 064 
u.uuz_ Q Q53 

-0 002+ ' 117 
u.uuz_ Q 144 


14.55l°; 3 2 2 


SNIa (SNLS) 


13+ 006 

13+ 013 
U - 1,3 -0.12 


14+0-03 
u - ±y: -0.03 

14+0-06 
u-1 -0.07 


20+ 16 
u - zu -0.35 


14.96li;^ 


SNIa(SNLS) 
+CMB+gas 


n i .7+0.01 

°- 17 -o.oi 

Q 17 +0.03 
O- 1 ' -0.04 


n i e + 0.01 

°- 16 -o.oi 

16+ ' 02 
u - ±D -0.02 


n nc + 0.05 
°- 05 -0.05 

0.051";- 


14 48+0-19 


SNIa(SNLS) 
+CMB+SDSS 
+gas 


22+ - 01 
u - z -o.oi 

n 99 +0.03 
u - zz -0.03 


n ic+0-01 

u - io -o.oi 
n 1 5+0-02 

u - ±o -0.02 


01+ ' 04 

u - ui -0.04 

Ooi+o.io 
u - u -o.io 


13.881°;^ 


SNIa(SNLS) 
+CMB+SDSS 
+LSS+gas 


21+ - 01 
u - zl -o.oi 

o 91+0 03 


n i c+0.01 

°- 16 -o.oi 

u - lo -0.02 


01+ ' 04 

u - ul -0.04 
01+°' 10 

u - ul -o.io 


13.881°;^ 



A dimensionless and independent of Hq version of SDSS observational parameter 



A = D v (z sdss ) 



^sdss 

-^0-Pl( 2 sds S ; Mm, j^rj 
H{z s dss'j ^m, ^r c )z^ dss (l + Z s d ss ) 2 



1/3 



(47) 
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where Dy (z s dss) is characteristic distance scale of the survey with mean redshift z s & ss 
70,73,74 . We use the robust constraint on the DGP model using the value of A = 
0.469 ± 0.017 from the Luminous Red Galaxy (LRG) observation at z s ^ ss = 0.35 [70] . 
This observation permits the addition of one more term in the \ 2 °f equations (I34|l 
and (|46[) to be minimized with respect to H(z) model parameters. This term is: 

2 _ [A>bs ~ Ah] 2 

XSDSS - 2 ( 4 °> 

sdss 

The baryon gas mass fraction for a range of redshifts is another observational 
test, can also be used to constrain cosmological models H(z). The basic assumption 
corresponding to this method is related to the baryon gas mass fraction in clusters |54| 
[75][76] as: 

this quantity is constant, related to the global fraction of the universe £2^ j f^rn • jSgas 
can be written as: 

q _ 1 M b _ b n b 

bsas - (l+/3)M tot - 0TpjJ2^ m 
where b is a bias factor suggesting that the baryon fraction in clusters is slightly lower 
than for the universe as a whole. Also 1 + /3 is a factor taking into account the fact that 
the total baryonic mass in clusters consists of both X-ray gas and optically luminous 
baryonic mass (stars), the later being proportional to the former with proportionality 
constant (3 ~ 0.19y/h [54] . Dimensionless parameter for this observation is given by 
[74]: 



i + p n m \D A (z-,n m ,Or ) 



^\D L (z;n m ,n rc )J { ' 

where D^ at is the angular diameter distance corresponding to flat pure CDM (f2 m = 1). 
Least square quantity in the likelihood analysis for this observation is: 

2 IC D D \ — V* [*^gas S ( z i) ~ Sg&s( z i \ ^m, f2r c , £)] 2 
X Iti "m, Mr a ) — / o , 

i 

(52) 

Marginalizing over the nuisance parameter, £ gives: 

x 2 sas (n m ,n rc ) = K-^- (53) 



where 



oobs/ \2 

gas ; J (54) 



W = 



qobs / 
-'gas ( 



( z i) ' ^gasfej ^m, f?r c , £ — 1) 



(55) 
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Fig. 8 Evolution of density contrast as a function of scale factor for different values of Q Tc in 
a flat universe. 



and 



Y = 2^ ~2 ( 56 ) 



We use the 26 cluster data for Sg&s( z ) reported in Ref. [54] to examine DGP modified 
gravity model. According to equations (J3U, (l4*Tj) . (gBJ, (gHJ and (JSHJ) we can con- 

strain free parameters of the model using observational data set related to background 
evolution. 

In what follows we perform a combined analysis of SNIa, CMB, gas cluster and 
SDSS to constrain the parameters of the DGP model by minimizing the combined 

2 2 2 2 2 

X = XSNIa + XcMB +Xgas + XsDSS- The best values of the model parameters from the 
fitting with the corresponding error bars from the likelihood function marginalizing over 
the Hubble parameter in the multidimensional parameter space are: Q m = 0.30^q'q 5 , 
Q rc = 0.12t°;° 2 and Q K = 0.01±°;^ at la confidence level with X m in/ N d.o.f = 0.94. 
The Hubble parameter corresponding to the minimum value of \ 2 ls h = 0.61. Here 
we obtain an age of 14.52^ ^ Gyr for the universe. Using the SNLS data, the best fit 
values of model parameters are: O m = 0.22^;^, Q Tc = 0.15±° ^ and Q K = O.Ol^g ^ 
at la confidence level with Xmin/Nd.o.f — 0.85. Table [2] indicates the best fit values 
for the cosmological parameters with one and two a level of confidence. 

Using the peaks position we find different values for the present matter density and 
Ta - Table [3] illustrates the best fit values and corresponding derived age of universe. 
According to the values reported in Tables [2] and [3] we infer that the value of f2 rc 
is very sensitive to the observational results from CMB. As usual we take the values 
confined using shift parameter of CMB instead of one given by absolute values from 
peaks position as a reliable results. 
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Table 3 The best values for the parameters of DGP self accelerating model with the corre- 
sponding age for the universe from fitting with SNIa from new Gold sample and SNLS data, 
SNIa+first peak and SNIa+first peak+SDSS+LSS also for third peak experiments at one and 
two a confidence level. 



Observation 




f2 rc 




Age (Gyr) 


SNIa(ncw Gold) 
+First Peak 


u - io -0.03 

0-15l o;°ol 


o.ostHI 


42+ - 10 
u -^ z -o.io 

n 42+ - 16 


15.86+°^ 


SNIa(ncw Gold) 
+First Pcak+ 
SDSS+LSS 


0-30±g;gi 

U - dU -0.04 


o-oil°;°4 

01+ ' 02 
u - u -o.oi 


n rry+0.08 

57+ 09 
u -°'_0.09 


1 3 77+0-09 
' -0.09 


SNIa(new Gold) 
+Third Peak 


13+ - 04 
u - io -o.oi 

1Q+0.08 
u - lo -0.02 


09+ 01 

U - UM _0.02 
09+ ' 02 


41+O.O8 
41+ - 16 




SNIa(ncw Gold) 
+Third Pcak+ 
SDSS+LSS 


o 2q+°- 01 

n 9q +0.02 


01+ ' 01 
u - ui -o.oi 

01+ ' 02 
u - u -o.oi 


o.58t° :S! 


13.851^9 


SNIa(SNLS) 
+First Peak 


n+ - 01 
u -- L1 -o.oi 

ii+°- 02 

u - 11 _0.02 


ii+ - 01 
u - ±1 -o.oi 

0-Hl 0°02 


n o fi +0.05 
Ucl0 -0.05 

n o fi +0.09 


±<i.»»_ 23 


SNIa(SNLS) 
+First Pcak+ 
SDSS+LSS 


u - io -o.oi 
0-15l°o:°03 


0.08+° °oi 

o.o8+° °2 


42+ - 05 
U - 4Z -0.05 

42+ 09 
u -** z _0.09 




SNIa(SNLS) 
+Third Peak 


o.o8l°:S 


13+ - 01 

n is+°' 02 

u - ±J -0.02 


33+ 04 

n oq+0.09 
u -*->_0.09 


15 ai+0.29 
lo - 91 -0 .31 


SNIa(SNLS) 
+Third Pcak+ 
SDSS+LSS 


16+ - 03 

u - lo -0.03 
u - iD -0.04 


07+ ' 02 
u - ul -0.01 

n n7+ - 03 


45+ - 10 

u -^°-0.09 

n 45+ - 14 


14 04+°' 40 
l^.U^i_ 43 



4 Constraints by Large Scale structure 

So far we have only considered observational results related to the background evolu- 
tion. In this section using the linear approximation of structure formation we obtain 
the growth index of structures and compare it with the result of observations by the 
2-degree Field Galaxy Redshift Survey (2dFGRS). 
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Fig. 9 Growth index versus redshift for different values of Q Ta . Here we imagined flat universe. 



Koyama and Maartens [44] have recently shown the evolution of density perturba- 
tions requires an analysis of the 5-dimensional gravitational field. In this model Poisson 
equation is modified and shows the suppression of growth due to gravity leakage. The 
continuity and modified Poisson equations for the density contrast 5 = 5p/p in the 
cosmic fluid provide the evolution of density contrast in the linear approximation (i.e. 
<5< l)[H[77l[78] as: 



5 + 2-5 
a 



v s V 



4ttG 



0, 



where 



^l-Mfl+J,) 



(57) 
(58) 



the dot denotes the derivative with respect to time. Thus the growth rate receives an 
additional modification from the time variation of Newtons constant through a. 

The effect of dark energy in the evolution of the structures in this equation enters 
through its influence on the expansion rate. The validity of this linear Newtonian 
approach is restricted to perturbations on the sub-horizon scales but large enough 
where structure formation is still in the linear regime [44U77j[78j . For the perturbations 

rV2, 



larger than the Jeans length, Aj 
matter (CDM) reduces to: 



G 1 + 



! p, equation (|57[1 for cold dark 



5 + 2-5 
a 



(59) 



The equation for the evolution of density contrast can be rewritten in terms of the 
scale factor as: 
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Numerical solution of equation (|60|) in the FRW universe in the background of DGP 
model is shown in Figure [8] In the CDM model, the density contrast 5 grows linearly 
with the scale factor, while we have a deviation from the linearity as soon as universe 
enters to acceleration era. Increasing Q Tc leads a decreasing in the evolution of density 
contrast which is in agreement to the finding about the behavior acceleration parameter 
versus Q Ta (see Figure [TJ. 

In the linear perturbation theory, the peculiar velocity field v is determined by the 
density contrast [771179] as: 

v(x) = H Q f f S(y)-^^d 3 y, (61) 

±* J |*-y| 

where the growth index / is defined by: 

f=^, (62) 
a in a 

and it is proportional to the ratio of the second term of equation (|59|) (friction) to the 
third (Poisson) term. 

We use the evolution of the density contrast 8 to compute the growth index of 
structure /, which is an important quantity for the interpretation of peculiar velocities 
of galaxies, as discussed in 79,80] for the Newtonian and the relativistic regime of 
structure formation. Replacing the density contrast with the growth index in equation 
(|60[) results in the evolution of growth index as: 



d f _ 3H ft , 1 



d In a 2a 



l 



aH 2 



(63) 



Figure [9j shows the numerical solution of equation (|63[) in terms of redshift. An analytic 
formula for the present growth index in the flat ylCDM model has been given in ref. 
[45] as f(z = 0.0, Q m ) ~ Q m e , here we extend this formula for the universe governed 
by DGP modified gravity. The simplest form for fitting formula in the wide range of 
Qm and Q Ta is 

f(z = 0.0; Q m , Qr c ) Q? n 6 + (0.0159 + 0.0603J? m ) x 

exp f [1.0694 - 0.3867 In Q m ] 2 Q r 

-0.0542« m - 0.0100 (64) 

Figure [TU] shows growth index, f(z = 0.0; Q m = 0.30), as a function of Q Ta derived 
from numerical solution of equation (|63[) and illustrated by fitting formula (equation 
(64}). 

To use observational results implied to linear structure formation we rely to the 
observation of 220, 000 galaxies with the 2dFGRS experiment provides the numerical 
value of growth index |70| . By measurements of two-point correlation function, the 
2dFGRS team reported the redshift distortion parameter of e — //re — 0.49 ± 0.09 at 
2 = 0.15, where re is the bias parameter describing the difference in the distribution of 
galaxies and their masses. Verde et al. (2003) used the bispectrum of 2dFGRS galaxies 
[8il82] and obtained K verde = 1.04 ± 0.11 which gave / = 0.51 ± 0.10. Now we fit the 
growth index at the present time derived from the equation (|63[) with the observational 
value. This fitting gives a less constraint to the parameters of the model, so in order 
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to have a better confinement of the parameters, we combine this fitting with those of 
SNIa+CMB+SDSS which have been discussed in the previous section. We perform the 
least square fitting by minimizing x 2 = xlrna + Xcmb + Xgas + xIdss + Xlss> where 

2 {fobsjz = 0.15) - fthjz = 0.15; Q m , Orjf f . 

XLSS — 2 < bt> ) 

a fa bs 

The best fit values with the corresponding error bars for the model parameters by using 
new Gold sample data are: O m = 0.28^° °^, Q Ta = 0.13±g$ and Q K = -0.002+°;°^ 
at ler confidence level with Xmin/^d.o.f = 0.93. Using the SNLS supernova data, the 
best fit values for model parameters are: Q m = 0.21^'^, Q Ta — 0.16^q'q^ and Qk = 
0.01+^4 at lcr confidence level with Xmin/ N d.o.f = 0.84. The error bars have been 

— 2 

obtained through the likelihood functions (£oc e x ' ) marginalized over the nuisance 
parameter h [S3]. The best values reported in the Ref. [31] using SNIa+CMB+SDSS 
are: Q m = 0.270, fl Tc = 0.125 for Gold sample SNIa and for SNLS SNIa are: Q m = 
0.255, Q Ta = 0.130, while in Ref. [26, using SNIa+SDSS, Q m = 0.270^\l, Q Ta = 
0.216^0 ^3. We concluded that observational results from large scale structure given by 
2dfGRS in addition to including baryon gas mass fraction results, put weak constraints 
on the DGP model free parameters. 

The likelihood functions for the three cases of (i) fitting model with Supernova data, 
(ii) combined analysis with the three experiments of SNIa+CMB+gas and (iii) combin- 
ing all five experiments of SNIa+CMB+gas+SDSS+LSS are shown in Figures [121 and 
1131 The joint confidence contours in the (n m , fir c ) plane are also shown in Figures ITU 
and[l5]for Gold sample and combined observational results, respectively. Figures [TBI and 
[TTlshow the joint confidence interval for SNLS data and SNIa+CMB+gas+SDSS+LSS 
experiments. 
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5 Age of Universe 

The "age crisis" is one the main reasons of the acceleration phase of the universe. The 
problem is that the universe's age in the Cold Dark Matter (CDM) universe is less than 
the age of old stars in it. Studies on the old stars [84] suggest an age of 13^ Gyr for 
the universe. Richer et. al. [85] and Hasen et. al. [86] also proposed an age of 12.7 ±0.7 
Gyr, using the white dwarf cooling sequence method (for full review of the cosmic age 
see [5])- The age of universe integrated from the big bang up to now is given by: 



Figure [TT] shows the dependence of I/nio (Hubble parameter times the age of universe) 
on Qr c for a flat universe. Obviously increasing Q Tc results in a longer age for the 
universe. As shown in the lower panel of Figure 1111 Q Tc behaves as the same as dark 
energy, Qa, in the flat ylCDM model. Finally we do the consistency test, comparing 
the age of universe derived from this model with the age of old stars and Old High 
Redshift Galaxies (OHRG) in various redshifts. Table [2]shows that the age of universe 
from the combined analysis of SNIa+CMB+gas+SDSS+LSS is 14.55t° 22 G y r and 
13.88^q ' 15 for new Gold sample and SNLS data, respectively. These values are in 
agreement with the age of old stars |84] , Here we take three OHRG for comparison 
with the DGP model, namely the LBDS 53W091, a 3.5-Gyr old radio galaxy at z = 1.55 
[87] . the LBDS 53W069 a 4.0-Gyr old radio galaxy at z = 1.43 [88] and a quasar, APM 
08279 + 5255 at z = 3.91 with an age of t = 2.lt^Gyr [B5]. The latter has once 
again led to the "age crisis". An interesting point about this quasar is that it cannot 
be accommodated in the /1CDM model [90]. To quantify the age-consistency test we 
introduce the expression r as: 



where t(z) is the age of universe, obtained from the equation (|66[1 and t b s is an 
estimation for the age of old cosmological object. In order to have a compatible age for 
the universe we should have r > 1. Table [4] shows the value of r for three mentioned 
OHRG. We see that the parameters of DGP model from the combined observations 
don't provide a compatible age for the universe, compared to the age of old objects, 
while the SNLS data result in a longer age for the universe. Once again for the DGP 
model, APM 08279 + 5255 at z = 3.91 has a longer age than the universe but gives 
better results than most Quintessence and braneworld models 91 ,92 ,93 94 . 

6 Conclusion 

We studied a self accelerating cosmological model, DGP modified gravity. The effect of 
this model on the age of universe, the radial comoving distance, the comoving volume 
element and the variation of the apparent size of objects with the redshift (Alcock- 
Paczynski test) have been studied. The evolution of density contrast, S, as a function 
of scale factor for various values of Q Tc shows that increasing fi Tc suppresses the 




to 



(66) 



r = 



t(z; n m ,Or c ) _ t(z; n m , Mr c )H a 



(67) 
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Fig. 11 -ffo^O ( a g e of universe times the Hubble constant at the present time) as a function 
of Q rc in a flat universe. Increasing Q rc gives a longer age for the universe. Lower panel shows 
the same function versus Six m the flat ylCDM model and to = —1.0. 



growth of density contrast, which is in agreement with the behavior of acceleration 
parameter versus Sl Tc ■ We extrapolate the relation of the growth factor in terms of 
Qi to the present time and showed that the power-law term is the dominant term the 
DGP model. To constrain the parameters of model we fit our model with the new 
Gold sample and SNLS supernova data, CMB shift parameter, position of the first 
and third peaks of power spectrum of temperature fluctuations at the last scattering 
surface, the Cluster Baryon Gas Mass Fraction, location of baryonic acoustic oscillation 
peak observed by SDSS and large scale structure formation data by 2dFGRS. The 
best parameters obtained from fitting with the new Gold sample data are: h = 0.62, 
Q m = 0.28±jj;^, Qr c = 0.13ig;{£ and Q K = -0.002+°;o53 at la confidence level with 
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Table 4 The value of r for three high redshift objects, using the parameters of the model 
derived from fitting with the observations. 



Observation 


LBDS 
53W069 
2 = 1.43 


LBDS 
53W091 
2 = 1.55 


APM 
08279 + 5255 
2 = 3.91 


SNIa (new Gold) 


, 07+O.I2 
-0.12 


1 19 +0.10 

1 - ±z -o.io 


7 o+0'07 
u - la -0.07 


SNIa(new Gold)+CMB 
+SDSS+gas 


i ot+0.03 
1,OJ -0.08 


1 11+0.03 
-0.07 


n 00+O.O2 


SNIa(new Gold)+CMB 
+SDSS+LSS+gas 


i 07+O.O5 
-0.03 


1 ,,+0.04 
1 ' lz -0.03 


U - S4 -0.03 


SNIa (SNLS) 


1 r^o + 0.17 
l.oo_ 23 


1 26+ 014 
i - ZD -0.19 


1 00+°' 13 
l.UU_ 21 


SNIa(SNLS)+CMB 
+SDSS+gas 


1 36+ 003 
1 -°°-0.03 


1 19 +0.02 
1,1 -0.02 


n oc+0.02 
U - St, -0.02 


SNIa(SNLS)+CMB 
+SDSS+LSS+gas 


1 36+ 003 
1 -°°-0.03 


1 19 +0.02 
1,1 -0.02 


n oc+0.02 



X 2 m in/ N d.o.f = 0.93 and by using the SNLS data are: Q m = 0.22^;^, Q rc = 0.15±g g^ 
and Qk = O.OIIq'q^ at la confidence level with Xmin/^d.o.f — 0.84. Comparing our 
results to that of previous results |26j[31j showed that large scale structure observations 
from 2dFGRS experiment had weak effect on confining the acceptance intervals for the 
free parameters. The observational constraint just by using SNIa+CMB indicated that 
our universe is spatially open but combining these result with SDSS+gas+LSS showed 
that our universe in the DGP model is very good agreement with the spatially flat 
universe. In comparison between ylCDM and DGP in terms of xi, Table [5] shows that 
these two models result almost same values. 

We also performed the age test, comparing the age of old stars and old high redshift 
galaxies with the age derived from this model. From the best fit parameters of the model 
using new Gold sample and SNLS SNIa, we obtained an age of 14. 55^22 Gyr and 
13.88^q']J Gyr, respectively, for the universe which is in agreement with the age of old 
stars. We also chose two high redshift radio galaxies at z — 1.55 and z — 1.43 with a 
quasar at 2 = 3.91. The ages of the two first objects were consistent with the age of 
universe, i.e., they were younger than the universe while the third one was not. 
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